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A theory of additive Markov chains with long-range memory, proposed earlier in 
Phys. Rev. E 68, 06117 (2003), is developed and used to describe statistical proper- 
ties of long-range correlated systems. The convenient characteristics of such systems, 
a memory function, and its relation to the correlation properties of the systems are 
examined. Various methods for finding the memory function via the correlation 
function are proposed. The inverse problem (calculation of the correlation function 
by means of the prescribed memory function) is also solved. This is demonstrated 
for the analytically solvable model of the system with a step-wise memory function. 
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I. INTRODUCTION 



The problem of long-range correlated dynamic systems 
a long time in many areas of contemporary physics 111 u 



RCS) has been under study for 

n nnn * nnnni^^^Wn iology ^ 

|ll|, |pj, economics |a |pj, |bj, literature [lj, U3, LL2I, La UM > etc - U U- One of the ways 
to get a correct insight into the nature of correlations in a system consists in constructing 
a mathematical object (for example, a correlated sequence of symbols) possessing the same 
statistical properties as the initial system. There exist many algorithms for generating 
long-range correlated sequences: the inverse Fourier transformation (20 



2l\ . the expansion- 
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modification Li method 
correlated Levy walks 



22j, the Voss procedure of consequent random additions 



2J, the 



24j . etc. |20j. The use of the multi-step Markov chains is one of the 



most important among them because they offer a possibility to construct a random sequence 
with necessary correlated properties in the most natural way. This was demonstrated in 
Ref. where the concept of Markov chain with the step-wise memory function was 

introduced. The correlation properties of some dynamical systems (coarse-grained sequences 
of the Eukarya's DNA and dictionaries) can be well described by this model 25 1. 

A sequence of symbols in the Markov chain can be thought of as the sequence of states 
of certain particle, which participates in a correlated Brownian motion. Thus, every L-word 
(the portion of the length L in the sequence) can be considered as one of the realizations 
of the ensemble of correlated Brownian trajectories in the "time" interval L. This point 
gives an opportunity to use the statistical methods for examining the correlation properties 
of the dynamic systems. Another important reason for the study of Markov chains is its 
application to the various physical objects [2^, 21, 3- e -S-> to the Ising chains of spins. The 
problem of thermodynamics description of the Ising chains with long-range spin interaction 
is still unresolved even for the ID case. However, the association of such systems with the 
Markov chains can shed light on the non-extensive thermodynamics of the LRCS. 

In this paper, we ascertain the relation between the memory function of the additive 
Markov chains and the correlation properties of the systems under consideration. We exam- 
ine the simplest variant of the random sequences, dichotomic (binary) ones, although the 
presented theory can be applied to arbitrary additive Markov processes with finite or infinite 
number of states. 

The paper is organized as follows. In the first Section, we introduce the general relations 
for the Markov chains, derive an equation connecting the correlation and memory functions 
of additive Markov chains, and verify the robustness of our method by numerical simulations. 
The second part is devoted to the study of the correlation function for the Markov chain 
with the step-wise memory function. In Subsec. flHI Bj) we reveal a band structure of the 
correlation function and obtain its explicit expression. Subsec. (fill CJ) contains the results 
of asymptotic study of correlation function. 



3 



II. GENERAL PROPERTIES OF ADDITIVE MARKOV CHAINS 

A. Basic notions 

Let us consider a homogeneous binary sequence of symbols, = {0, 1}, i G Z = 
—2, —1, 0, 1, 2, .... To determine the N-step Markov chain we have to introduce the con- 
ditional probability P(cij | aj_jv, Ui-N+i, ■ ■ ■ > a i-i) of the definite symbol a, (for example, 
aj = 1 or a,j = ) occurring after the iV-word T^, where Tjy j denotes the sequence of 
symbols ai-N, a«-Af+i, • • • , flj-i. Thus, it is necessary to define 2^ values of the P-function 
corresponding to each possible configuration of the symbols in the iV-word T/v,i- Since we 
intend to deal with the sequences possessing the memory length of the order of 10 6 , we need 
to make some simplifications. We suppose that the P-function has the additive form, 

N 

P(a i = l\ T Ntl ) = Y,f(a i - k ,k). (1) 
fc=i 

Here the value f(ai-k, k) is the additive contribution of the symbol ai-k to the conditional 
probability of the symbol unity occurring at the zth site. Equation corresponds to the 
additive influence of the previous symbols on the generated one. Such Markov chain is 
referred to as additive Markov chain, Ref. j^]. The homogeneity of the Markov chain is 
provided by the independence of the conditional probability Eq. ([TJ of the index i. It is 
possible to consider Eq. as the first term in expansion of conditional probability in the 
formal series of terms that correspond to the additive (or unary), binary, ternary, and so on 
functions up to iV-ary one. 

Let us rewrite Eq. in an equivalent form, 

N 

P{ ai = 1 | T Nti ) = a + ^P(r)(a,_ r - a). (2) 

r=l 

Here 

TV N 

d = J2mr)/[l-Y,(f^r)-f(0,r))] 

r=l r=l I I 

is the average number of unities in the sequence, Ref. |29j, and 

F(r) = /(l,r)-/(0,r). 

We refer to P(r) as the memory function (MF). It describes the strength of impact of 
previous symbol aj_. r upon a generated one, a,. Evidently, this function has to satisfy 
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condition ^ P(cti = 1 | T N ^) ^ 1. To the best of our knowledge, the concept of the 
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25 1 where it 



memory function for multi-step Markov chains was introduced in papers 
is shown that it is convenient to use it in describing the correlated properties of complex 
dynamical systems with long-range correlations. 

The function P(a,i = 1 | Ijv,i) contains complete information about correlation properties 
of the Markov chain. In general, the correlation function and other moments are employed 
as the input characteristics for the description of the correlated random systems. Yet, the 
correlation function takes account of not only the direct interconnection of the elements cij 
and cii +r , but also their indirect interaction via other elements. Our approach operates with 
the "origin" characteristics of the system, specifically with the memory function. 

The positive values of the MF result in persistent diffusion where previous displacements 
of the Brownian particle in some direction provoke its consequent displacement in the same 
direction. The negative values of the MF correspond to the antipersistent diffusion where 
the changes in the direction of motion are more probable. In terms of the Ising long-range 
particles interaction model, which could be naturally associated with the Markov chains, 
the positive values of the MF correspond to the attraction of particles and negative ones 
conform to the repulsion. 

We consider the distribution Wi{k) of the words of definite length L by the number k of 

L 

unities in them, ki(L) = Y2 a i+h an d the variance D(L) of k{(L), 

i=i 



D{L) = {k-k)\ (3) 



where the definition of average value of g(k) is g(k) = ^2 g(k)WL(k). 

k=0 

Another statistical characteristics of random sequences is the correlation function, 

K(r) = a { a i+r - a 2 . (4) 

By definition, the correlation function is even, K(—r) = K(r), and K(0) = a(l — a) is the 
variance of the random variable <2j. The correlation function is connected with the above 
mentioned variance by the equation 

K(r) = ^(D{r-l)-2D(r)+D(r + l)), (5) 



or 



in the continuous limit. 



ld 2 D(r) 
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B. Derivation of main equation 

In this subsection we obtain a very important relation connacting the memory and corre- 
lation functions of the additive Markov chain. Let us introduce the function 0(r) = p(a« = 
1 | aj_ r = 1), which is the probability of symbol <jj = 1 occurring under condition that 
previous symbol aj_ r is also equal to unity. This function is obviously connected to the 
correlation function K(r), see Eq. (@J) since the quantity aiai- r is the probability of simul- 
taneous equality to unity of both symbols, and aj_ r . It can be expressed in terms of the 
conditional probability <p(r), 

a-a— = acj)(r). (7) 
Substituting Eq. into Eq. (J3J we get: 

K(r) = a<f)(r) - a 2 . (8) 

For the TV-step Markov chain, the probability of the symbol aj = 1 occurring depends 
only on the previous iV-word. Therefore, to obtain the value of 0(r) one needs to average 
the conditional probability P Eq. (J2J) over all realizations of the iV-words with the condition 

Q'i—r 1; 

0(r) = p(a,i = 1 | a,i- r = 1) (9) 
= ^P(aj = 1 | T N)i )p(T N>i | di- r = 1). 

If the value of r is less or equal to N, then ai- r in Eq. (jHJ) is one of the symbols 
Oj_i, aj_2, . . . , a,i-N hi the word Tjv,i- In this case, the sum in Eq. © allows not all iV- 
words, but only the words that contain the symbol unity at the (i — r)th position. If r > N, 
the memory function F(r) equals zero in this region and, therefore, the sum in Eq. © 
contains all terms corresponding to all different iV-words. 
Substituting Eq. Q into Eq. Q, we obtain: 

^(r) = o^P(T jy|i |o j _ r = l) (10) 

N 

+ 2j F(r') ^(ai_ r / - a)P(T N:i \ a,i- r = 1). 

r'=l Tjv.i 
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According to the normalization condition, the first sum in Eq. (fTUj) is equal to unity. 
Consider the sum 

22 a i-r'P{ T N,i I 0-i-r = 1) (11) 

in the second term in RHS of Eq. (JTOj) . The symbol aj_ r / is contained within the word Tjv,i- 
Therefore, Eq. (fTT|) represents the average value of aj_ r / under condition a^ r = 1. In other 
words, it equals the probability <ft(r — r') of a$_ r / = 1 occurring under the condition aj_ r = 1: 

ai- r >P(T N ,i | a,_ r = 1) = <f>(r - r'). (12) 

w 

Substituting this equation into Eq. ([TD|) we obtain: 

JV 

+ J^F(r')(0(r-r')-a). (13) 



iV 

f /• i ■ f/ 

r'=l 



Taking into account Eq. (|SJ), we arrive at the relation between the memory function and the 
correlation function: 

N 

K{r) = ^F(r')K(r-r'), r ^ 1. (14) 

□ 

This equation was first derived by the variation method in Ref. [29J. 

Another equation resulting from Eq. (JT4"Jl by double summation over index r establishes 
a relationship between the memory function F(r) and the variance D(L), 

N 

M{r, 0) = F(r')M(r, r'), (15) 

r'=l 

M(r, r') = D(r - r') - (D(-r') + r[D(-r' + 1) - D(-r')}). 

Equation (jSJ) and parity of the function D(r) are used here. 

The last equation shows convenience of using the variance D(L) instead of the correlation 
function K(r). The function K(r), being a second derivative of D(r) in continuous approx- 
imation, is less robust in computer simulations. It is the main reason why we prefer to 
use Eq. (I15jl for the long-range memory sequences. This is our tool for finding the memory 
function F(r) of a sequence using the variance D(L). 
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C. Numerical reconstruction of the memory function 

Let us verify the robustness of our method by numerical simulations. We consider a 
model memory function, 

'l-r/10, l^r<10, 



F(r) = 0.1 



(16) 



0, r > 10. 

shown in the inset in Fig. ^ by solid line. Using Eq. we construct a random unbiased 



0.03 



0.02 



0.01 



o.oo 




Figure 1: Calculated correlation function K(r) of the Markov chain constructed with the model 
memory function F(r), Eq. I|16|l. shown by the solid line in the inset. Dots in the inset correspond 
to the memory function reconstructed by solving Eq. (|14[) with correlation function K (r) presented 
in the main panel. 

a = 1/2, Markov chain. Then, with the help of the constructed binary sequence of the 
length 3 x 10 6 , we calculate numerically the correlation function K(r) by solving the set of 
N linear equations ([14)1 . The result of these calculations is given in Fig.0 One can see that 
the correlation function K(r) mimics roughly the memory function F(r) over the region 
1 ^ r ^ 10. In the region r > 10, the memory function is equal to zero but the correlation 
function does not vanish J^. Then, using the obtained correlation function K(r), we solve 
numerically Eq. (j!4)l . The result is shown in the inset in Fig. Qby dots. One can see an 
excellent agreement of initial, Eq. ()16|) . and reconstructed memory functions F(r). 

The ability of constructing a binary sequence with an arbitrary prescribed correlation 
function by means of Eq. (JUJ) is the very nontrivial result of this paper. 

Yet another approach to numerical finding the memory function is an iteration procedure. 
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For its realization let us rewrite Eq. ([14)1 in the form, 



F(r) = ^- f ^^V). (17) 
{ ) K(0) ^ K(0) { ' K ' 

Using Eq. (|17|) with starting iteration F (r) = 0, we obtain the formula, 

, s ff(r) A K(r - r') , A 
WO = ^ " S^Jf.M, » > < 18 > 

Thus, the memory function can be presented as the series, 

(j *(o) £ ^ 2 (o) + r ^ r ^ k 3 (o) tiyj 

Note, that the Markov chain with the definite correlation function K(r) exists if the se- 
ries f|19|) is convergent and the obtained function implies the probability Eq. (J2J) satisfy- 
ing the requirement ^ P{fli = 1 | T^,i) ^ 1 for arbitrary word Tn,i- If a = 1/2 we 

obtain the restriction, El-^( r )l ^ 1- The sufficient, but not necessary, requirement is 
N-i ~^ 
£ \K(r)\^ 1/12 -\K(N)\/3. 

r=l 



III. CORRELATION FUNCTION OF THE CHAIN WITH THE STEP- WISE 

MEMORY FUNCTION 

In the previous section, we obtained the relationship (|14j) between two characteristics of 
the Markov chain, the memory and correlation functions, and used this equation to solve the 
problem of finding the memory function via the known correlation function. Here we present 
the solution of the inverse problem. We suppose the memory function to be known and find 
the correlation function of the corespondent Markov chain. To simplify our consideration, 
we examine the step-wise memory function, 

(a, r < N, 

F(r) = (20) 
I 0, r> N. 

The restriction imposed on the parameter a can be obtained from Eq. (J2J): \a\ < 1/N. Note 
that each of the symbols unity in the preceding N-word promotes the emergence of new 
symbol unity if < a < 1/N. This corresponds to the persistent diffusion. The region of 
parameter a, determined by inequality — 1/N < a < 0, corresponds to the anti-persistent 
diffusion. If a = 0, one arrives at the case of the non-correlated Brownian motion. 
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A. Main equation for the correlation function 

Substituting Eq. (|2Uj) into Eq. (jT4*)l we obtain the relation, 

TV 

K{r) = a^K(r -r'), r ^ 1. (21) 

r'=l 

Here, the correlation function is assumed to be even, K{—r) = K(r). Equation (}2"Tj) is 
the linear recurrence of the order of N for r ^ N + 1, so we stand in need of N initial 
conditions. For the unbiased sequence, a = 1/2, we have K(0) = 1/4. The solution of 
Eqs. (J21|) written for r = 1, . . . , N yields the constant value of the correlation function, 
K(r) =K atr = 1,...,N- 1, 

k> - 4(i -„(*-!» • (22) 

Subtracting Eq. (|2 lj) from the same equation written for r + 1, we derive another, more 
convenient, form of the recurrence: 

K(r + 1) - (1 + a) K(r) + a K(r - N) = 0. (23) 

This equation is of the order of N + 1, thus we need an additional initial condition. It can 
be derived from Eq. (}2Tj) : K(N) = Kq. Note that the possibility to rewrite Eq. (}2Tj) in the 
form of Eq. ()23|) is the result of the simple structure of the memory function. We solve the 
obtained recurrent equations by the most natural method, by means of step-by-step finding 
the sequent values of the correlation function. Such an approach is very suitable for the 
analysis of correlation function at r > N. 



B. Correlation function at r > N 

1. Band structure of the correlation function 

Equation (J2"T|) allows one to find numerically the unknown correlation function K(r). The 
result of this step-by-step calculation is presented in Fig. by solid line. One can easily see 
the discontinuity of K(r) at the point L = N = 100, the breakpoint of the curve is observed 
at L — 2N. Such behavior of the correlation function is the result of using the memory 
function of the step-wise form. To clarify this fact it is convenient to change the variable r 
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Figure 2: Correlation function K (r) (solid line) obtained by two different methods: the numerical 
simulation of Eq. (|14|) and exact solution Eq. 1)411) . The dotted line is for the contribution to K(r) 



of the first root of Eq. (|32|) . The dashed line refers to the correlation function obtained in 
The vertical lines indicate the limits of the bands numbered by s. 



IS]- 



by the band number s and the intra-band number p: 

K(r) = K s (p), r = sN + p+l, p = 0, 1, . . . , N - 1, s = 0, 1, . . . . (24) 

Within sth band, Eq. (f2"H|) is the recurrence of the second order with the term a K s -i(p) that 
is determined at the previous step, while finding the correlation function for the (s — l)th 
band. 



2. General expression for the correlation function 

In the zeroth band (s = 0, 1 ^ r ^ N), as it was shown above, the correlation function 
is constant, 

K (p) = K . (25) 

For the first band (s = l, N +1 ^ r ^ 2N), taking into account that K(r — iV — 1) = Ko(p), 
we have 

K^p) = K (1 - (1 - aN)(l + a)"). (26) 

The correlation function decreases quasi- continuously within the first band. However, as it 
was mentioned above, there exists a discontinuity in the K(r) dependence at r = N. This 
discontinuity disappears in the limiting case of the strong persistence, a — > 1/JV. 
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Substituting the obtained formula (|2^|) in Eq. (j2HJ), we find the solution K 2 (p) for the 
second band (s = 2, 2iV + 1 ^ r < 3iV), 

tf 2 (p) = K (1 - (1 - aN)((l + a) p+N - pa(l + af' 1 ). (27) 

The correlation function K(r) is continuous at the interface between the first and second 
bands, Ki(N) = fG(0). However, its first derivative of K{r) is discontinuous here (see 
Fig. |2J). Using the induction method, one can easily derive the formula for K s (p) in the sth 
band (sN + 1 < r < (s + 1)N): 

K s {p) =K Q U-{\- aN) EC-a) 4- ^! + a)^ N+ ^ +1 C^ )N+ ^ , (28) 

C k rjn + i) 

n r(jfe + i)r(n - k + 1) ' 

It follows from Eq. ([28)1 . that the first (s — 2) derivatives of the correlation function K(r) 
are continuous at the border between the (s — l)th and sth bands, but the derivative of the 
(s — l)th order changes discontinuously. Under the condition aN <C 1, Eq. (|28|) takes a 
simpler form, 

X s (p) = Jjfoo-C^.^^. (29) 

It is seen that the correlation function decreases proportionally to a s with an increase of the 
band number s. 

It is not easy to analyze the asymptotical behavior of the function K(r) at large s because 
the number of summands in Eq. (|2*Hj) increases being proportional to s. It is the reason to 
propose another approach for the asymptotical study of the correlation function K(r) at 
s > 1. 



C. Asymptotical study of the correlation function 

1. Derivation of the characteristic equation 

The general solution of linear recursion equations (|21|) can be represented as the linear 
combination of iV different exponential functions, 

N 

K{r) = Y< a &- ( 30 ) 
i=i 



12 



To find the values of £j, we substitute the fundamental solution, 

K(r) = f , (31) 

into Eq. (J2T|) and obtain the characteristic polynomial equation of the order of N. Constant 
multipliers ai are to be determined by initial conditions. 

It is more convenient to use Eq. (|2*3*|) instead of Eq. (J2T|) . that implies the characteristic 
equation of the order of N + 1, 

£N+1 _ (1 + a ^N + a = Q (32) 

The extra root of this equation, £ = 1, appears as a consequence of passing on to the 
equation of order of iV + 1 from that of the order of N . The corresponding coefficient, eij, 
in Eq. ()3()j) is equal to zero because the correlation function should decrease at r —>■ oo. 

Our study shows that Eq. (|3*2|l has one real positive root less than unity in the case of 
odd N. In the case of even N, there are two real roots, one positive and one negative. The 
remaining roots are complex. All absolute values of roots are less than unity, that is in 
agreement with the finiteness of memory function F(r). In the case of large N, the absolute 
magnitudes of all roots are close to unity for nearly all values of a satisfying the inequality, 

llni«l. (33) 

Distribution of the roots in the complex plane £ is shown in Fig. El 

1,0 

0,5- 

0,0- 
-0,5- 
- 1,0- 

- 1,0 - 0,5 0,0 0,5 1,0 

Figure 3: The dots are the roots of the characteristic equation Q32|) for N = 100 and a = 0.008. 
The solid line is the circle |£| = 1. 




In the simplest case, iV = 



2, Eq. has two real roots, 
a [a? 

6, 2 = 2 ± VT + a ' 



(34) 
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Taking into account the initial conditions, we find the solution of Eq. ()21j) in the form, 

cv 

K{r) = — -_=(E-i(l - fc) - - 6))- (35) 

4(1 — aJV cr + 4a 

This expression can be simplified at small and large values of parameter a. For a C 1, one 
obtains 

K{r) = -^ r ^ +1 (36) 

with square brackets standing for the integer part. The correlation function in the sequent 
odd and even points are equal to each other. In accordance with Eq. (J29|) . K(r) decreases 
at r — > oo being proportional to a s . In the opposite limiting case of the strong persistency, 
a — > 1/2, we have two different roots, 

6 = 1-^0, 6 = -^ + ^, (37) 

with = 1/2 — a. The coefficient corresponding to the second root is much less than that 
corresponding to the first one. Besides, the second term in Eq. (|35|) decreases more rapidly. 
Therefore, the approximate solution in this case is 

K(r) = - exp -40(r - l)/3. (38) 

2. Correlation function at small a 

Let us return to the case of arbitrary value of N. If a is very small, i.e. at 

ilnl»l, (39) 

Eq. has iV roots with small absolute magnitudes: 

& = a 1/A Vs(27r^) + i sm(2vr ^)), = 0, . . . , JV - 1. (40) 

The correlation function, being a linear combination of the power functions with these roots 
as their exponents, decreases proportionally to a s , which agrees with Eq. (|2~5|l. 

The coefficients tij in the linear combination Eq. (J3*U|) can be found in a general case, 
without any restrictions imposed on the value of a. The solution of Eq. (J21j) written for 
1 ^ r ^ N — 1 along with K(0) = 1/4 can be expressed with the help of the Vandermond 
determinants: 

K(r) = K (aN - 1) £ , (41) 

k=i n (e fc -e,-) 

with = 1. 
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3. Correlation junction at not too small a 

In the case flSlIj) of not too small a, the absolute magnitudes of all roots are close to unity. 
It is convenient to rewrite Eq. (jH2J), introducing two new real variables 7 and tp instead of 
complex x according to 

X = (1 - -L)e*. (42) 



Equation takes the form, 



a(eH^-l) = (l--L)e^-l. (43) 
For the real root, Eq. yields, 

aiV 7 (e 1/7 - 1) = 1. (44) 

This expression along with Eq. (|31|) determines the asymptotical behavior of correlation 
function. It was first obtained in Ref. Q|. The qualitative approach of this paper did not 
allow one to take into account the contribution of other roots as it is done in the present 
paper. 

Equation (|4~3*|) yields all remaining complex roots with the values of (p, which are quite 
uniformly distributed over the circle [0, 2tt] and 

-~ln-. (45) 
7 a 

The roots of Eq. ()32j) located in the vicinity of point £ = 1 are shown in Fig. 0] The single 
real root is much closer to the line Re £ = 1 than the other ones. Besides, the coefficients 
Oj in Eq. (jjlj) (see also Eq. (|3()|l) for all terms containing the complex exponents are much 
less than those for the term with the real exponent. Therefore, the behavior of correlation 
function K(r) is determined generally by the term with the real exponent. 

The exact correlation function K(r) resulting from the numerical simulation of Eq. (|41|) 
and its approximation determined by the contribution of the real root alone are shown in 
Fig. 121 These curves are compared with that obtained in Ref. by a qualitative method. 

The obtained correlation function can be used to calculate one of the most important 
characteristics of the random binary sequences, the variance of number of unities in the 
L-word. The results of the numerical simulations are shown in Fig. One can see a good 
agreement of curves plotted using both of these methods. 
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Figure 4: The roots of characteristic polynomial equation close to the point £ = 1 for iV = 4000, 
a = 2 • 10 -4 . The solid line is Re £ = 1. 
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Figure 5: Variance D{L) for the Markov chain with N = 100, a = 0.008 obtained by means of 
exact Eqs. (|41|) . ([5]) (solid line) and using only one root of Eq. (|32[) (dashed line). The thin solid 
line describes the non-correlated Brownian diffusion, D(L) = L/4. 

4- Conclusion 

Thus, we have demonstrated the efficiency of describing the symbolic sequences with 
long-range correlations in terms of the many-step Markov chains with the additive memory 
function. Actually, the memory function appears to be a suitable informative "visiting card" 
of any symbolic stochastic process. Various methods for finding the memory function via the 
correlation function of the system are proposed. Our preliminary consideration shows the 
possibility to generalize our concept of the Markov chains on larger class of random processes 
where the random variable can take on arbitrary, finite or infinite number of values. 

The suggested approach can be used for the analysis of different correlated systems in 
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various fields of science. For example, the application of the Markov sequences to the theory 
of spin chains with long-range interaction makes it possible to estimate some thermodynamic 
characteristics of these non-extensive systems. 
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